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ABSTRACT: This paper presents a modified PSO algorithm that utilizes the PSO with double chaotic 
maps to solve the Economic Ordered Quantity (EOQ) problem with different cost functions. In proposed 
approach, we employ the PSO method that involves the alternating use of chaotic maps in estimating the 
velocity of the particle. Presently the PSO method has been widely used in complex designs problems with 
multiple variables where the optimization of a cost function is required. The use of pure analytical method 
for the Economic Ordered Quantity (EOQ) problem is widely studies however the best solution may 
require problem dependent derivations, complex mathematical calculations and at the worst it fails to 
achieve global best solution, hence to overcome the problem this paper presents a different approach by 
using PSO with irregular velocity updates which is performed by some kind of random function which 
force the particles to search greater space for best global solution. However the random function itself 
derived from a well-defined mathematical expression which limits its redundancy hence in the paper we 
are utilizing the two different chaotic maps which are used alternatively this mathematically increased the 
randomness of the function. The simulation of the algorithm verifies the effectiveness and superiority of 
the algorithm over standard algorithms. 
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I. INTRODUCTION 

The deterioration of goods in present inventory system is an unavoidable problem, which is applicable 
from common goods like, foods, vegetables and fruit to the specialized goods like medicines. Hence it is 
important to consider the deterioration loss during the EOQ problem. Although the a number of deteriorating 
inventory models have been studied and proposed in recent years but the consideration of cost change due to 
ordering size has not been studied in best of our knowledge. Hence this paper presents a modeling and solution 
approach for such systems. Furthermore these proposal tackles the EOQ as the optimization problem and solved 
using numerical and analytical solution approach like lambdaiteration method, quadratic programming and the 
gradient method. In order to use any of these methods for EOQ should be capable of handling the 
constrains(cost change due to ordering size) which are inherently highly non-linear this leads to multiple local 
minimum points of the cost function. In order to overcome these limitations this paper presented the PSO 
technique enhanced by alternative use of two different chaotic maps. The rest of the paper is arranged as second 
section presents a brief review of the related literatures in third section discusses the problem formulation while 
fourth section explains the PSO and the variants used in the paper finally in chapter five and six respectively 
presents the simulated results and conclusion. 



II. LITERATURE REVIEW 

In recent years, a lots of research have been performed on EOQ models and evolutionary systems some 
of them found useful during writing this paper are presented here. Lianxia Zhao [7] studied an inventory model 
with trapezoidal type demand rate and partially backlogging for Weibull-distributed deterioration items and 
derived an optimal inventory replenishment policy. Kuo-Lung Hou et al. [10] presents an inventory model for 



I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 5 I Iss. 8 I August 2015 I 69 I 



A Modified PSO Based Solution Approach for Economic Ordered Quantity Problem with Deteriorating 



deteriorating items considering the stock-dependent selling rate under inflation and time value of money over a 
finite planning horizon. Liang YuhOuyanget al. [1] presented an EOQ inventory mathematical model for 
deteriorating items with exponentially decreasing demand. Their model also handles the shortages and variable 
rate partial back ordering which dependents on the waiting time for the next replenishment .Kai- Wayne Chuang 
et al. [2] studied pricing strategies in marketing, with objective to find the optimal inventory and pricing 
strategies for maximizing the net present value of total profit over the infinite horizon. The studied two variants 
of models: one without considering shortage, and the other with shortage .Jonas C.P. Yu [4] developed a 
deteriorating inventory system with only one supplier and one buyer. The system considers the collaboration 
and trade credit between supplier and buyer. The objective is to maximize the total profit of the whole system 
when shortage is completely backordered. The literature also discuss the negotiation mechanism between 
supplier and buyer in case of shortages and payment delay. The model allows shortages and partially 
backlogging at exponential rate. Ching-Fang Lee et al. Al [14], considered system dynamics to propose a new 
order system for integrated inventory model of supply chain for deteriorating items among a supplier, a producer 
and a buyer. The system covers the dynamics of complex relation due to time evolution. Michal Pluhacek et al 
[15] compared the performance of two popular evolutionary computational techniques (particle swarm 
optimization and differential evolution) is compared in the task of batch reactor geometry optimization. 
Both algorithms are enhanced with chaotic pseudo-random number generator (CPRNG) based on Lozi chaotic 
map. The application of Chaos Embedded Particle Swarm Optimization for PID Parameter Tuning is presented 
in [16]. Magnus Erik et al [17] gives a list of good choices of parameters for various optimization scenarios 
which should help the practitioner achieve better results with little effort . 

III. PROBLEM FORMULATION 

The objective of an EOQ problem is to minimize the total inventory holding costs and ordering costs 
which should be found within the equality and inequality constraints (operational constrains) limitations. The 
simplified cost function of each inventory item can be represented as described in (2) 

n 

C T = £<:,($) (3.1) 

i = 1 

q(S;) = * S t (3.2) 

where 

C T = Total Inventory Cost 
q = Cost Function of Invertory i 
a t = per unit cost of Invertory i 
S t = Ordered size of Invertory i 

3. 1.1 Order Size Dependent Cost Model 1: The inventory with order dependent cost shows a greater variation in 
the total inventory holding cost and ordering cost function. Since the cost depends upon order size, a cost 
function contains higher order nonlinearity .The equation (3) presents a cost function quadratic behavior: 

C i (S i ) = a i +l3 i S i +Y i S? 3.3) 



where 



a it fii and y t are the cost coefficient of inventory i, 

And can be explained as 

a is the minimum cost and should be a positive constant 
ft is the price per unit and should be a positive constant 
y is the discount and should be defined as follows 
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Let a single inventory system 

C(S) = a + pS + yS 2 (3.4) 

Since the cost must always be higher for higher quantities even after adding discount hence it should be a 
monotonic increasing function in the ranged E [S min S max ]. Because the minimum ordered quantity must be 
equal to or larger than 1 we can say that the 

dC(S) r ~i . 

— > 0, 5 E [S m [ n S max ] (3.5) 

P + 2 yS >0, S E [S min S max ] 

B 

Y > ~2S* Se l S rnin Smax ] (3-6) 

The equation (6) states the condition required for monotonic function however in the present case on more 
consideration is required because the cost per unit should be decreased when ordered at higher amount but the 
total cost should never fall below the total cost of lowered ordered quantity as shown in (7). 

C(5j) > C (S]), for i > j, and {5), 5}} e [S min S max ] ... (3.7) 



Or 



d 2 C(S ) r . , N 

^2 < 8 [ S minSmax ] ( 3 - 8 ) 

2y<0, Se[S min S max ] 
y< o, se[s min s max ] (3.9) 



Combining the equation (6) and (9) we can get the possible values for the y as follows 

~27s <r<0 ( 3 - 10 ) 




s 

y= 0.0S y= 0.06 y= 0.05 y= 0.04 



Figure 1: the effect of y on cost C(5), where a = 5, /? = 10, and S max = 100. 
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3.1.20rder Size Dependent Cost Model 2: this model represents the model where items per unit cost depends 
upon order size segment which can be represented with several piece wise quadratic functions reflecting the 
effects of different inventory ordering costs according to order size segment it belongs to .In general, a 
piecewise quadratic function is used to represent the input-output curve of an inventory with lumped order size 
[4] and described as 

Let the inventory is divided into N equal segments then the size of each segment. 

S — S’ 

c max u min 

S seg ~ ~ 

Now the quantities under the j th segment 

Sj = {Smin + O' - 1) * S Seg >S min + 0 - 1) * S seg + L 

Smin + 0 - 1) * S seg + 2 , S min + 0) * S seg ~ 1} (3.11) 

Pj Be the cost per unit when ordered under j th segment size, hence the ordering cost 

Cj=Pj*Sj (3.12) 



Th e/5y should be estimated as 



Pi < Pj , for i > j 



And 



C t > Cj,for i > j 

Now the maximum number of quantity that can ordered under the j th segment 

S = S -+■ 7 s|e S — 1 

‘-’j ,max u min 1 J u seg ± 

While the minimum quantity that can ordered under the (j + l) th segment 

S(j +1 ),min S m i n + (0 "f 1) 1) * S se g 

Now the cost for Sj >max and S(j +1 ^ min must follow 

S(j +l),min * P(j+1) ^ Sj.max * Pj >Pj+ 1 < Pj 
(S m in + (O' + 1) - l) * S seg ) * P( J+ i) > (S min + j* S seg - 1) * Pj 



Smin * Pj+1 + j * S seg * Pj +1 > S min * Pj + j * S seg * Pj - Pj 
Smin (Pj+1 — Pj) + j * ^eg (Pj+1 — Pj) > — Pj 
(Pj+1 - Pj)(S min + j * S seg ) > -Pj0 



Pj+1 

Pi 



< 1 



Pj+1 <Pj 1 - 



(^min "f j * ^seg ) 
1 



(^min "f j * S se g ) 



.(3.13) 
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Figure 3: shows the cost per unit for different segments (/?&)( above) and cost of ordered quantity under 
different segments for = 10,5^^ = 5 ,S max = 95 ,S seg = 6, N = 15. 

In general, total inventory ordering costs is calculated by equation (4). Determining the selection of total 
inventory order quantity for each unit is dictated by the demand and backlogging rate, and can be solved by 
economic total inventory holding and ordering cost. This paper assumes that such selection is given a-priori. 
Therefore, to obtain an accurate and practical EOQ solution, the total inventory holding and ordering cost 
function should be considered with both multi-item total inventory holding and ordering costs with order size 
segment effects simultaneously [7]. Thus, the total inventory holding and ordering cost function (3) should be 
combined with (4), and can be represented as follows: 
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Where 



Ct(.S t ) = 



CnVd, 

CaiSi), 



segment 1, 
segment 2, 



S’. . 

Lj i,min 

Si, i 



< Si < S U1 

<Si< s U2 






c ik (Si) = a ik + b ik Si + c ik sf (3.13) 



and S i min is the minimum order size of nventory i. 
St, max i s th e maximum order size of nventory i 



(3.12) 



3.2. Equality and Inequality Constraints 

3.2.1 Demand and Stock Balance Equation: For Demand and Stock balance ,an equality constraint should be 
satisfied . The total stock should be equal or greater than the total demand plus the total Deterioration loss 

n 

Si'demand. *^i,/oss (3.14) 

i = 1 

where S i demand S i)ioss represents the total demand and 
Deterioration loss of i th inventory is a function of the 
units ordered that can be represented using Demand (DJ 
and Deterioration (LQ coefficients [2] as follows : 

n n 

i = 1 i = 1 



Y S t L t (3.15) 




IV. MATHEMATICAL MODELING 

The mathematical model in this paper is rendered from reference [1] which derived the following conclusive 
equations: 

Notation Used: 

c 1 : Holding cost, ($/per unit)/per unit time. 
c 2 : Cost of the inventory item, $/per unit. 
c 3 : Ordering cost of inventory, $/per order. 
c 4 : Shortage cost, ($/per unit)/per unit time. 
c 5 Opportunity cost due to lost sales, $/per unit. 
t x Time at which shortages start. 

T Length of each ordering cycle. 

W The maximum inventory level for each ordering cycle. 

S The maximum amount of demand backlogged for each ordering cycle. 

Q The order quantity for each ordering cycle. 

Inv(t ): The inventory level at time t. 



The inventory level at any given time (t) ,t 1 < t <T: 



lnv{f) = y (ln[l + 8{T - t)] - ln^l + 8(T - U)]} (4.1) 

o 



Maximum amount of demand backlogged per cycle as follows: 



D 

S = - Inv(T ) = — ln[l + 8(T - tj] 
o 



.(4.2) 
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The ordered quantity per cycle is given by 



Q = W + S = -f— [e (e - A)t i - 

6 — A 



D 



l]+-ln[l + <J(7’-t 1 )] (4.3) 



The inventory holding cost per cycle is 

c x A 



HC = 



6(6 -X) 
The deterioration cost per cycle is 



0 -Xt x 



Q -I 

e 0tl - 1 - ~ Ol (4.4) 



DC = c 2 A 

The shortage cost per cycle is 

SC = c a D | 



- 1) - j(l - ef l£ )} ( 4 - 5 ) 



T - t x 1 



5 8 2 

The opportunity cost due to lost sales per cycle is 



ln[l + S(T - ti)]} ....(4.6) 



BC = c s D {(7 - - ^ln[l + S(T - tj]} (4.7) 



Considering the scenario 3.1.1 the cost is given by the equation considering single inventory type 

C = a + /? * S + y * S 2 (4.9) 

Now replacing c 2 in equation (4.5) by (4.9) 

DC = a + (3 * A j- — — — l) — - (l — e _/itl )| + y 

* A {trrj (e (e_A)tl - 1) - J (1 - e~ Xt f (4.10) 

Considering the scenario 3.1.2 the cost is given by the equation considering single inventory type 

Cj=Pj*Sj (4.11) 

Now replacing c 2 in equation (4.10) by (4.11) 

DC = - !) - - e-^fforAyfj(e (e ~^ - l) - j(l - e"^)} 

£ Sj (4.12) 

Therefore, the average total cost per unit time per cycle can be expressed as: 

TVC = TVC(t lf T ) 

= (holding cost +deterioration cost + ordering cost +shortage cost+ opportunity cost due to lost sales)/ length of 
ordering cycle 

For scenario 3.1.1 
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—At i 



TVC = 



1 f c x Ae 
T\6(6 - X) 



,6>ti 



1 e(e^-l) 



+ c 3 + C 4 (D) 



+ c 5(^) ( T ~ ti — 



+ yA 



ln(l + 5(7 - tO) 



+ cr + f3A 



7 - ^ in(l + 5(7 - tO) 
52 

,(e-A)ti _ ]_ 1 _ , 



, — Ati 



0-A 



^ j _ g-foi 



e -a 



.(4.13) 



For scenario 3.1.2 



TVC 



I c x Ae~ 



T \ 6(6 - X) 



Atl ( 6(e Xtl — l)\ /7 — Zq ln(l + 6(T — tf)) 

e Bt! _ X 1 ') + C 3 +C 4 (D)' 1 V 1JJ 



s 2 



( Zn(l + 5(7 — t0)\ ( e ( e 2^1 — l 1 — e Xtl 

+ c 5 (D)[T-t 1 ^ £ W- + pA 



6 - X 



. . (4.14) 



The objective of the model is to determine the optimal values of Zqand Tin order to minimize the average total 
cost per unit time, TVC. The optimal solutions t { and 7* need to satisfy the following equations: 

For scenario 3.1.1 



d TVC _ 1 
dt x _ 7l 



CiAAe Xtl f Qu Q{e Xtl — l)\ c x Ae Xt ^(Qe etl — 6e Xtl ) 

' 0(0 - A) V 1_1 X ) + 0(0 - A) 

+ cm (-5- 1 + 5(l + (J | r _ ti)) ) + CsW (-1 + (1 + 5(7 - tO)" 1 ) 



+ PA(t 



= 0 . 



(O-A)t 1 _ p -At\ 



) + 2yA 



f 6 -X)t 1 _ y y 



-A tl 



6 -A 



^e-A) tl _ e - 






. (4.15) 



and 



d TVC 1. . , 

^=f( c ‘ (D H s *(i + 



5(7 - tl ))) + C5(Z)) “ (1 + 5(7 ' tl))_1 )] 

1 (c x Ae~ Ul ( au 1 0(e Al 






7 2 \ 0(0 - A) 



1) 



+ c 3 + c 4 (D) 



7 - t 4 in(l + 5(7 - tO) 



5 2 



. Zn(l + 5(7 — tO)\ fe^ e 1 - )tl — 1 1 — e 

+ c 5 (D) 7-ti * * — ) + a + pA' 



-Ati 



+ yA 



>(8—x)t i _ j ^ 



6 — 2 . 



= 0. 



6 - X 

(4.16) 
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For scenario 3.1.2 



d TVC _ 1 



CiAXe Xtl ( 0ti i 0(e Atl - l)\ qAe Atl (0e etl - 0e Atl ) 



0(0 -A) 

+ C4 (D) ( —5 1 + 



A 

1 



6(1 + 5(7 - td) 

+ /?>l(e (e - A)t i - e“ X£l ) 1 = 0 (4.17) 



0(0 - A) 

+ c 5 (D) (-1 + (l + 5(7 — ti)) ) 



dTVC 1 . , . , 

57 _ r( C4 ^( 6 5(1 + 5(7 



^ _ ti)) ) + c 5 (0) (l - (1 + 6(7 - ti))' 1 ) j 



1 ( c x i 4e Atl / 0(e Atl — l)\ /7 — t x ln(l + 6(7 — tf)) 

' 1 ' e 0tl — 1 — ; I + c 3 + c 4 (D) 1 ^ 






7 2 1 0(0 - A) 

, / ln(l + 6(7 - tO) , 

+ c 5 (D) 7 - j — + PA 



5 6 2 

e (0 — A)ti _ y y _ e -Xti> 



0 — A 



= 0. 



.(4.18) 



V. PARTICLE SWARM OPTIMIZATION (PSO) 

The PSO algorithms inspired by the natural swarm behavior of birds and fish. In PSO each particle in the 
population represents a possible solution of the optimization problem, which is defined by its cost function. In 
each iteration, a new location (combination of cost function parameters) of the particle is calculated based on its 
previous location and velocity vector (velocity vector contains particle velocity for each dimension of the 
problem). 

Initially, the PSO algorithm chooses candidate solutions randomly within the search space. The PSO algorithm 
simply uses the objective function to evaluate its candidate solutions, and operates upon the resultant fitness 
values. 

Each particle maintains its position, composed of the candidate solution and its evaluated fitness, and its 
velocity. Additionally, it remembers the best fitness value it has achieved thus far during the operation of the 
algorithm, referred to as the individual best fitness, and the candidate solution that achieved this fitness referred 
to as the individual best position or individual best candidate solution. Finally, the PSO algorithm maintains the 
best fitness value achieved among all particles in the swarm, called the global best fitness ,and the candidate 
solution that achieved this fitness, called the global best position or global best candidate solution. 

The PSO algorithm consists of just three steps, which are repeated until some stopping condition is met: 

1 . Evaluate the fitness of each particle 

2. Update individual and global best fitness’s and positions 

3. Update velocity and position of each particle 

4. Repeat the whole process till the 

The mathematical equations involved in above operations are as follows: 

The velocity of each particle in the swarm is updated using the following equation: 

v(f + 1) = w * v(f) + c 1 * ( pBest — x(0) + c 2 * (, gBest — x(Q) (5.1) 

Modified PSO with chaos driven pseudorandom number perturbation 
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v(i + 1) = w * v(i) + c 1 * Rand * ( pBest — x(0) + c 2 * Rand * (, gBest — x(Q) (5-2) 

A chaos driven pseudorandom number perturbation (Rand) is used in the main PSO formula (Eq. (13)) that 
determines new “velocity’ ’and thus the position of each particle in the next iterations (or migration cycle). The 
perturbation facilities the better search in the available search space hence provides much better results. 

Where: 

v(i + 1) — New velocity of a particle. 
v(i) — Current velocity of a particle. 
c lf c 2 — Priority factors. 
pBest — Best solution found by a particle. 
gBest — Best solution found in a population. 

Rand — Random number, interval (0, 1). Chaos number 
generator is applied only here. 



x(i) — Current position of a particle. 

The new position of a particle is then given by (5.3), where x(i + 1) is the new position: 

x(i + 1) = x(i) + v(i + 1) (5.3) 



Inertia weight modification PSO strategy has two control parameters w start and w end . A new w for each 
iteration is given by (5.4),whereistandforcurrentiteration number and n for the total number of iterations. 



((.Wstart ~ W e nd ) * 0 
W = w start 



.(5.4) 



VI. CHAOTIC MAPS 

A chaotic map is an evolution function) that shows some sort of disordered behavior. These maps may 
be discrete-time or continuous-time. The terms is came from chaos theory where the behavior of dynamical 
systems which greatly depends upon the initial conditions. In such cases even a small change in initial 
conditions may cause the exceptionally diverging results which makes the long-term prediction for such system 
impossible .Presently there are many such maps have been discovered however in our work we have selected 
the Lozi and Logistic maps because of their simplicity. The logistic map is a polynomial mapping (equivalently, 
recurrence relation) of degree 2, shows the complex, chaotic behavior from very simple non-linear dynamical 
equations. Mathematically, the logistic map is written as: 

X n+ i=liX n (l-X n ) (6.1) 

The Lozi map is a simple discrete two-dimensional chaotic map presented by following equations: 

X n+1 = l-x\X n \+bY n (6.2a) 

Y n+ i=X n (6-2 b) 



VII. IMPLEMENTATION OF IMPROVED PSO ALGORITHMFOR ECONOMIC 
ORDERED QUANTITY (EOQ) PROBLEMS 

Since the decision variables in EOQ problems are t x and T with S = [S lr S 2 , ... . ,S n } where S t ordering 
quantity of i th inventory, the structure of a particle is composed of a set of elements corresponding to 
th e[t lf T,S]. Therefore, particle’s position at iteration k can be represented as the 

vector Xf = (P^, P^, , P^) where m = n + 2 and n is the number of inventories. The velocity of particle 

^corresponds to the generation updates forall inventories. The process of the proposed PSO algorithm can be 
summarized as in the following steps. 
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1. Initialize the position and velocity of a population at random while satisfying the constraints. 

2. Update the velocity of particles. 

3. Modify the position of particles to satisfy the constraints, if necessary. 

4. Generate the trial vector through operations presented in section 4. 

5. Update and Go to Step 2 until the stopping criteria is satisfied. 




Figure 8: Flow Chart of the Proposed Algorithm. 



VII. SIMULATION RESULTS 

The proposed IPSO approach is applied to three different inventory systems explained in section 3 and 
evaluated by all three PSO models as follows: 

• The conventional PSO 

• The PSO with chaotic sequences 

• The PSO with alternative chaotic operation 

The simulation of all algorithms is performed using MATLAB . The population size N P and maximum 
iteration number iter max are set as 100 andlOO, respectively. w max and w min are set to 0.9 and 0.1 respectively 
because these values are widely accepted and verified in solving various optimization problems. The list of all 
values used for the system are shown in the table below 
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Table 1: parameter values used for different PSO algorithms 



Name of Variable 


Value Assigned 


Cl 


2 


c 2 


1 


W max 


0.9 


W min 


0.1 


[i (logistic map) 


4.0 


k (logistic map) 


0.63 


a (lozi map) 


1.7 


b (lozi map) 


0.5 


Total Particles 


100 


Maximum Iterations 


100 



Table 2: values of system variables: 





Variable Value 


Variable Name 


Scenario 1 


Scenario 2 


A 


12 


12 


0 


0.08 


0.08 


s 


2 


2 


X 


0.03 


0.03 


c± 


0.5 


0.5 


c 2 


1.5 


1.5 


C 3 


10 


10 


c 4 


2.5 


2.5 


c s 


2 


2 


D 


8 


8 


a 


5 


5 


p 


10 


10 


r 


0.04 


N/A 


** mm 


1 


1 


^ max 


100 


95 


N 


N/A 


15 


^ seg 


N/A 


6 


Pi 


N/A 


10 




Figure 9 (a): surface plot for the normal inventory system. With respect to t x (the time at which 
shortage starts) and ^(ordering cycle time) the figure shows a smooth and continuous curve and hence can be 
solved by analytical technique also .Figure 9(b): surface plot for the order dependent inventory cost type model. 
With respect to t x (the time at which shortage starts) and ^(ordering cycle time) the figure shows much steep 
variations and continuous curve and hence can be solved by analytical technique with difficulties .Figure (c): 
surface plot for the order segment dependent inventory cost type model. With respect to fq (the time at which 
shortage starts) and T (ordering cycle time) the figure shows much abrupt variations and many discontinuities in 
the curve and hence can be very difficult to solve by analytical techniques. 
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Figure 10 (a): the value of objective function (fitness value or TVC) at every iteration of PSO for model 1. 
Figure 10(b): the best values of variables t 1 and T for all three PSO for model 1. 




Figure 11(a): the value of objective function (fitness value or TVC) at every iteration of PSO for model 2. 
Figure 1 1(b): the best values of variables t x and T for all three PSO for model 2. 




Figure 12: the value of objective function (fitness value or TVC) at every iteration of PSO for model 3. Figure 
12(b): the best values of variables t x and T for all three PSO for model 3. 



Table 3: Best Fitness Values by all three PSO for model 1, 2 and 3. 



Type of PSO 


Best Fitness (TVC) 


Model 1 


Model 2 


Model 3 


PSO 


11.6625 


17.7682 


14.9668 


PSOl 


11.4125 


17.5094 


14.7168 


PS02 


11.2736 


17.3695 


14.5779 



VIII. CONCLUSION AND FUTURE SCOPE 

In this paper presents the mathematical model for three different inventories systems where the 
inventory cost depends on order size by some continuous function of discontinuous function the paper also 
presents the derivations for evaluation of the function parameters for practical applications and finally it 
proposes an efficient approach for solving EOQ problem under the mentioned constrains applied 
simultaneously. Which may not be solved by analytical approach hence the meta-heuristic approach has been 
accepted in the form of standard PSO furthermore the performance of standard PSO is also enhanced by 
alternative use of two different chaotic maps for velocity updating finally it is applied to the EOQ problem for 
the inventory models discussed above and tested for different systems and objectives. The simulation results 
shows the proposed approach finds the solution very quickly with much lesser mathematical complexity. The 
simulation also verifies the superiority of proposed PSO over the standard PSO algorithm and supports the idea 
that switching between different chaotic pseudorandom number generators for updating the velocity of particles 
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in the PSO algorithm improves its performance and the optimization process. The results for different 
experiments are collected with different settings and results compared with other methods which shows that the 
proposed algorithm improves the results by considerable margin. 
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